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Abstract 


This paper describes several methods for the prediction of jet noise. All but one of 
the noise prediction schemes are based on Lighthill’s or Lilley’s acoustic analogy while 
the other is the jet noise generation model recently proposed by Tam and Auriault . 1 
In all the approaches some assumptions must be made concerning the statistical prop- 
erties of the turbulent sources. In each case the characteristic scales of the turbulence 
are obtained from a solution of the Reynolds-averaged Navier Stokes equation using a 
k — e turbulence model. It is shown that, for the same level of empiricism, Tam and Au- 
riault’s model yields better agreement with experimental noise measurements than the 
acoustic analogy. It is then shown that this result is not because of some fundamental 
flaw in the acoustic analogy approach: but, is associated with the assumptions made 
in the approximation of the turbulent source statistics. If consistent assumptions are 
made, both the acoustic analogy and Tam and Auriault’s model yield identical noise 
predictions. The paper concludes with a proposal for an acoustic analogy that provides 
a clearer identification of the equivalent source mechanisms and a discussion of noise 
prediction issues that remain to be resolved. 


1 Introduction 

The prediction of jet noise has been the object of continuous interest and study since the 
introduction of the jet engine for commercial use. The earliest theoretical formulation for 
aerodynamic noise was the work of Lighthill 2,3 . Lighthill’s equation was the first example 
of an “acoustic analogy.” The definition of an acoustic analogy, to be used in this paper, 
is any aerodynamic noise theory in which the equations of motion for a compressible fluid 
are rearranged in a way that seprates linear acoustic propagation effects. By definition, this 
rearrangement results in a set of equivalent sources that are assumed to be non-negligible in 
a limited region of space. In an apparent departure from formulations based on the acoustic- 
analogy, Tam and Auriault 1 have recently developed a jet noise prediction method in which 
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the sound sources are modeled explicitly and the propagation of sound from these sources 
is described by solutions to the linearized Euler equations. These two aerodynamic noise 
theories appear to be very different. In this paper we reconcile these apparent differences 
and show that if consistent assumptions are made concerning the statistical properties of the 
turbulent noise sources, both approaches can yield identical noise predictions. It should be 
noted that this reconciliation has only been demonstrated at 90 degrees to the jet axis. The 
different approaches could lead to very different results at other angles. 

In order to make predictions of the radiated noise it is necessary to describe the properties 
of the turbulence. Since a complete simulation of the jet turbulence and the noise it generates 
and radiates are too computationally expensive for high Reynolds number jet flows, noise 
predictions are often based on the solution of the Reynolds-averaged Navier Stokes (RANS) 
equations using a k — e turbulence model. Such a solution provides an estimate of the 
amplitude of the turbulent velocity fluctuations as well as a local length scale. Since the k — e 
solution only provides time-averaged properties, it is necessary to make assumptions about 
the statistical characteristics of the turbulence in order to predict the noise radiation that is 
an inherently unsteady phenomenon. In particular, the axial two-point cross correlation of 
the turbulent sources must be approximated. This correlation may be described in either a 
fixed frame of reference or a reference frame moving with the turbulence. 

One of the earliest attempts to couple estimates of the statistical properties of the turbu- 
lence from a steady flow prediction with a noise model based on the acoustic analogy was the 
MGB approach developed by Mani, Gliebe and Balsa: see Balsa and Gliebe. 4 More recent 
extensions have used RANS solutions for the flow field based on a k — e turbulence model. 
This has been referred to as the MGBK method: see Khavaran et al. 5 and Khavaran. 6 In 
developing a solution to the acoustic analogy equations these methods assume a form for the 
two-point cross correlation function for the turbulent sources in a moving frame of reference. 
Conversely, Tam and Auriault, 1 who also use RANS solutions for the flow field based on a 
k — e turbulence model, describe the two-point cross correlation in a fixed reference frame. 
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It is shown in this paper that, for the same RANS solution, jet noise predictions made with 
the Tam and Auriault model 1 give much better agreement with experimental measurements 
at 90 degreees to the jet axis than methods based on the aeosutic analogy. However, it will 
also be shown that this is not due to any inherent flaw in methods based on the acoustic 
analogy: but, is associated with the assumptions made concerning the statistical properties 
of the turbulent sources. Both approaches yield identical noise predictions at 90 degrees to 
the jet axis if consistent descriptions of the turbulent sources are chosen. 

This paper is organized as follows. First a noise prediction formula based on Lighthill’s 
acoustic analogy is developed. It is denoted here as Model I based on the acoustic analogy. 
The features of Tam and Auriault’s noise prediction model 1 are then described. Another 
model based on the solution to Lighthill’s equation is then formulated using the solution 
procedure of Tam and Auriault. 1 This is designated as Model II based on the acoustic 
analogy. The reasons for the differences between the noise predictions obtained with the 
different approaches are discussed and these differences are reconciled. Finally, an alternative 
form of acoustic analogy is proposed. It is argued that this form allows for the easier 
identification of the equivalent noise source mechanisms. 

2 Lighthill’s Acoustic Analogy: Model I 

In this section we develop a prediction scheme based on the solution to Lighthill’s equation. 
The details of the methodology do not follow those of the MGBK approach exactly. This is 
because the focus of the present paper is on noise radiation at 90 degrees to the jet axis, where 
mean flow/acoustic interaction effects are negligible. It is at this angle that the specification 
of the source and its assumed relationship to the k — e solutions is best assessed. Also, it 
is intended to keep the assumptions made and the level of empiricism used as consistent as 
possible between the different schemes described in this paper. So, the possible importance of 
the effects of anisotropy of the turbulence on the radiated noise, as proposed by Khavaran, 6 
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is not included. The analysis in this section follows that given by Goldstein' and Lilley in 
Chapter 4 of Hubbard. 8 It is repeated in sufficient detail here to emphasize any assumptions 
made and the differences with the alternative approaches presented in later sections of the 
paper. 

LighthiH’s equation 2 may be written in Cartesian tensor form as, 

d 2 p' 2 dV _ m 

dt 2 0 dxid.i'i dxjdxj 

where Ty is the Liglitliill stress tensor given by 

Tij = pUiUj + Sij [{p - p 0 ) - c 2 0 {p- p 0 )\ (2) 

and primes denote perturbations about the basic state denoted by a subscript o. c 0 is a 
constant speed of sound that is sensibly taken to be the speed of sound in the uniform 
medium surrounding the source region. ?/, is the instantaneous velocity vector. Viscous 
terms have been neglected in the Lighthill stress tensor. In the subsequent analysis it is 
assumed that the departures from isentropic behavior are everywhere small and that the 
How is at relatively low Mach number. Since we are concentrating on noise radiation at 90° 
to the jet axis it is also assumed that the primary contributions to the Lighthill stress tensor 
involve products of velocity fluctuations. The terms that are linear in the fluctuations on 
the right hand side of Eqn. (1) should be regarded as terms associated with the propagation 
of the sound and be placed on the left hand side of the equation: see Lillev. 9 These effects 
are negligible at 90° to the jet axis so this assumption is reasonable here. So we approximate 
the Lighthill stress tensor by 


T 7 / / 

ij — PsUplj 


( 3 ) 


where p s is the mean density in the source region. In the far field the density fluctuation is 
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readily shown to be given by 


(x ' () = i^y/ // m T ™ { yA - j V 1 ) dy (4) 

v(y) 

where T xx is the component of the Liglithill stress tensor in the direction of the far field 
observer and x = |x — y| & |x|. 

The far field spectral density for the intensity is related to the Fourier transform of the 
autocorrelation function of the far field pressure. 


S (x, u) 


OO 

J_ f {P f (^ t )'P f (^ t + T )) c iWT dr 
2 tt J PoC-o 

— oo 


(5) 


where ( ) denotes an ensemble average. Since, in the far field, p r = c 2 0 p ' we obtain 


S( *’“ )= 32 Adix> j .1 J (tF (6) 

-oo V(y, i V(y 2 ) 

where 


1 1 

h 


t — 


x-yi 

Co 


t + T 


x - y 2 


(7) 

( 8 ) 


If the turbulent statistics are assumed to be stationary and the usual far field approximation 
is made we obtain, 

OO 

s ^ ^ = 3 ^; C 5 X 2 J J J ( 1 ™ (yi^) (y2, T>)) e'“ T dyic/y 2 dr (9) 

-oo y(yi) y(y 2 ) 
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where 


To 


t+T + *. fo-y'> 

x c 0 


( 10 ) 


The two-point cross correlation function of the Lighthill stress tensor in a fixed reference 
frame may be denoted by, 


= (T xx (yi, t) T xx (y 2 , t + r)) 


(11) 


where rj = y 2 - yi, then, 


S (x, w) = 


LU 


32tt 3 p 0 c^x 2 


R f (y 1 ,r/,r)exp iw 


-°° V'(yi) V(t)) 


T — 


x 77 


X Co 


dy 1 dr] dr (12) 


As noted by Lighthill 2,3 and others, it is best to include as many properties of the source 
as possible prior to any modeling of the turbulent sources. To include the effects of source 
convection the statistical properties of the sources may be described in a moving frame of 
reference. This also has the advantage that it is the temporal variation in this frame that 
controls the noise radiation. In a fixed reference frame, the temporal variation is dominated 
by convection effects. For example, as noted by Goldstein,' a frozen pattern of turbulence 
conveeting subsonically would radiate no noise. However, its local time variation would 
depend on the convection velocity and the turbulent length scales and would not be zero. 
Let 


£ = ri- i c 0 M c t 


( 13 ) 


where, i is a unit vector in the direction of the mean flow and M c is the convection Mach 
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number of the turbulent eddies. This gives 


S (x, lj ) = 


cu 


327 i z p 0 c b o x 2 

-oo v( yi )v(a 
Rm(yi,£,r)d£dyidT 


exp < iuj 


(1 — M c cos 9)r — 


x c n 


x 


( 14 ) 


where R m denotes the two-point cross correlation function of the Lighthill stress tensor in 
the moving reference frame. Also, 


cos 6 = .r i fx 


(15) 


The wavenumber /frequency spectrum of the turbulent sources is given by 

OO 

H{y 1 ,a,u) = — ^-4 J J e t( ‘ 1JT ~ a<) R rn {y 1 ^ 1 T)cl^clr (16) 

V(€) -oo 

where a is a wavenumber vector. This describes the spatial and temporal periodicity of the 
source. Then, 


S (x, u) 


7 TUJ 4 1 

2 p 0 c 5 0 x 2 


[ H 

yi, — ,w (1 - M c cos9) 

1 

xc a 


1 (y> ) 


dyi 


(17) 


This shows that the far field noise depends on the components of the source wavenum- 
ber/frequency spectrum with a wavenumber that gives a sonic velocity in the direction of a 
far field observer and at a Doppler shifted frequency. 

To this point, other than the far field assumption, no approximations have been made. 
However, to proceed further, it is necessary to introduce a model for the two point cross 
correlation. It is usually assumed that, in the moving frame of reference, the correlation 
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takes on a Gaussian form (see Ffowcs Williams 10 and Lilley 8 ), 


Rm(yi,£,r) 


A 2 p 2 s u A s exp 



2 2 
-UJ S T 


( 18 ) 


where, i s is a characteristic length scale, u s is a characteristic frequency in the moving frame, 
u s is a velocity scale that characterizes the turbulent velocity fluctuations, and A determines 
the magnitude of the correlation. Here, it assumed that the characteristic length scale is the 
same in all directions. This restriction could easily be relaxed. It should be noted that this 
is simply a model for the turbulent statistics in a moving reference frame. It is not an exact 
relationship. Then, 


H \ yi, — ,w (1 - M c cos 6 ) 
1 xc„ 


A2 


(2tt) 


4 Ps U s 


UJc 


exp 


oj 2 (i — m c cos oy 
4 u 2 


exp 


(«Gf 


where, a is the magnitude of the wavenumber vector. Now, 


Oit s 


COs^-s 

C 0 


u s 

~ r^j 777 , 

C 0 


(19) 


( 20 ) 


where m is a characteristic Mach number for the turbulence and provides a measure of the 
compactness of the source region. For compact sources, m « 1, so that at 90 degrees to 
the jet axis we obtain, 

S( * w, = 12^/ ***•*(£) exp {“£f} rfyi (21) 

v 

If a RANS k — e solution is available it is possible to determine the contribution of each 
elemental volume in the numerical grid to the radiated noise spectrum. Here, k and e are 
the turbulent kinetic energy and visous dissipation rate per unit mass respectively. From 
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equation (21) this contribution is given by 


dS (x, oo) = 


A 2 


327 rp 0 c b 0 x 


.5^2 r ! 8 ! s 


00 


P»U s t' s Ws — ex P - 


00 c 


00 

4uj 2 


dV 


( 22 ) 


The length and time scales may be obtained from the k — s solution. We assume that 


u> s = 2 tt/t,, r s = c 7 (k/e), and 0 = ce{k^ 2 /f). 


(23) 


The k — e solutions indicate that, along the location of maximum shear, both r s and i s vary 
nearly linearly with axial distance. Then, with u s = y/2k/3, 


dS (x, co) = 


A 2 d 


7 T 


4 9 p 0 c b o x 2 


p 2 ,k 7,2 


00 


00c 


exp 


00 

4^2 


dV 


(24) 


In this form there are only two combinations of constants that may be determined by com- 
parison with experiments. The factor 


oo 

— I exp 
CO s 


00 

4u s 2 


(25) 


describes how each volume element contributes to a range of frequencies about the local 
characteristic frequency oo s . 

Figure (1) shows a comparison of the predicted radiated noise with experimental data 
by Tanna et al. 11 The one-third octave experimental data have been converted to spectral 
density assuming a smooth spectrum. The jet is operating at M ] = 0.911 and Tj/T 0 = 0.975. 
The jet diameter is 0.0508 m and the observer location is at 90° to the jet axis at a distance 
of 72 jet diameters. The k — e solution has been obtained using the code developed by Thies 
and Tam. 12 The grid is described by Thies and Tam. 12 It grows in physical size in a stepwise 
manner as the solution is marched in the axial direction. For the present calculations the 
solutions are saved at every quarter of a jet diameter downstream. The volume of the 
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elements is determined by the radial grid spacing, with a constant axial spacing and an 
axisymmetric assumption. The coefficients used in the noise predition are c T = 1.43 and 
A 2 e 3 = 1.93 x 10 -2 . The predicted spectrum has been matched to the measured spectrum 
at the peak frequency. The predicted frequency variation of the spectrum increases as / 2 at 
low frequency and decreases as f ~ 2 at high frequencies. The much faster decay at very high 
frequencies is due to the lack of resolution in the flow prediction near the nozzle exit. It 
is clear that the behavior of the experiments is quite different. The experimental variation 
of the spectral density increases approximately with / at low frequency and decreases as 
/ -1 at high frequencies. Figure (2) shows how each axial slice of length Dj of the 

jet contributes to the total spectrum. The spectral shape of the contribution of each axial 
slice is determined by Eqn. (25). The slices near the jet exit provide the high frequency 
contributions to the spectrum and the peak frequencies of the spectra decrease monotonically 
in the axial drection. It is clear that this spectral distribution is unable to provide sufficiently 
high levels away from the peak, particularly at low frequencies. 

3 Tam and Auriault’s Method 

In a recent development, Tam and Auriault 1 developed a model for noise radiation from 
“small-scale” turbulence. As the name suggests, this model, they argued, accounts for the 
noise radiation by the small-scale components of the turbulence. In an earlier paper, Tam 
et al. 13 had shown how experimental noise data for a wide range of jet operating conditions 
could be correlated extremely well with two shapes of spectral density function. One, with 
a well-defined peak and relatively rapid decay at high and low frequencies, was associated 
by Tam et al. 13 with noise generated by the large-scale structures in the turbulence. This 
spectrum matched the experimental noise data at small angles to the jet axis (including the 
peak noise direction). The second spectrum shape was much broader, similar to the measured 
spectrum shown in Fig. (1). This they associated with noise from the small scale turbulence 
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and it matched the measured spectrum at larger angles to the jet axis. Tam and Auriault 1 
argue that their noise prediction model provides a description of the latter mechanism. It 
should be noted that there is no direct evidence that these two separate mechanisms are 
actually responsible for the total noise radiation: but, the experimental correlations and the 
success of the model can not be overlooked. 

Before providing some of the details of Tam and Auriault’s analysis, it is useful to summa- 
rize the differences between the model of the previous section based on the acoustic analogy 
and Tam and Auriault’s model as well as the differences in the development of a prediction 
formula for the far field noise. Firstly, the Lighthill equation is based on the full equations 
of motion whereas Tam and Auriault’s model contains a heuristic argument to describe the 
noise source and uses the linearized Euler equations to describe the propagation of sound 
generated by the model sources. In the model based on the acoustic analogy it is assumed 
that the source is compact and the two-point cross correlation function for the source is 
modeled in a moving reference frame. The Green’s function in the model based on the 
acoustic analogy is simply the free space Green’s function for the wave equation. In Tam 
and Auriault’s method 1 there is no explicit assumption concerning the compactness of the 
source, (though compactness is implied at one point), the two-point cross correlation func- 
tion is formed in a fixed frame of reference, and the Green’s function is obtained from the 
adjoint solution of the linearized Euler equations. A final difference between the two formu- 
lations is that the far field noise depends on a model for the cross correlation of the Lighthill 
stress tensor in the model based on the acoustic analogy and on the cross correlation of the 
convective derivative of the source term in Tam and Auriault’s model. The consequence of 
this difference is examined in a later section of this paper. 

In Tam and Auriault’s model 1 it is proposed that the small scale turbulence generates a 
local pressure fluctuation that is proportional to the local turbulent kinetic energy per unit 
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volume. Tam and Auriault note that, in the kinetic theory of gases, the pressure is given by 


p = —mn (v • v) = ~p(v 2 ) (26) 

o o 

where v is the random molecular velocity, p is the density of the gas, and ( ) denotes an 
ensemble average. Invoking a direct analogy, the small scale turbulence is considered as 
small blobs of fluid that interact in each element of the turbulent flow. They argue that the 
resulting pressure is related to the kinetic energy of the fine-scale turbulence. That is, 

1 2 

Pturb Qs gP (v ) g Mb (27) 


where k s is the kinetic energy of the small-scale turbulence per unit mass. It is argued that 
once sound is generated by these local pressure fluctuations, the propagation of the noise 
may be described by the linearized Euler equations with a “source” term on the right hand 
related to the fluctuating pressure gradient generated by the fine-scale turbulence. That is, 


P 


d u[ _ du[ 


, dui 
u,- 


3 dx ; 


dp' 

dxj 


dq s 

dxi 


(28) 


where, an overbar denotes a mean quantity and u[ and p' denote acoustic field variables. The 
sound waves are also assumed to satisfy the linearized equations of energy and continuity 
and the equation of state for a perfect gas. That is, 


dp 

dt 


U.i 


dp 

dx i 


- dui 

7Px— = 0 

dx) 


(29) 


Tam and Auriault 1 then write the component equations in cylindrical polar coordinates 
(r, </>, x) and assume that the mean flow is parallel. 

Rather than converting the operators on the left hand sides of Eqns. (28) to (29) into the 
Lillev equation, 9 Tam and Auriault 1 make use of the adjoint equation. This is a very elegant 
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approach to describe the mean flow acoustic interaction effects. A complete development for 
three-dimensional mean flows is given by Tarn and Auriault. 14 Tam and Auriault 1 show that 
the periodic Green’s function for the linearized Euler equations is related to the solution of 
the adjoint Euler equations by, 

Pi(x 0 ,xi,w) = u a (x x , x 0 , CO ) 

p 2 (x 0 ,x i,w) = /•„ (X|.X„.>') (30) 

p3 (x 0 , X x , LO ) = W a (xi , X 0 , UJ ) 

p ti for n l. 2, 3 are the Green’s functions for sources in the x, r, and 0 components of 
the linearized momentum equations respectively and. u a , v a , and w a are the solutions to the 
adjoint linearized momentum equations. Note the reciprocal dependence on x 0 and X| . The 
adjoint equat ions may be written, 


p < lLUV a 


_dp, 

I'jJpa + U 


f)x 


p < tuu a 


u 


-dug 

dx 


dv n clu 


dx dr 


u ■ 


P < OjJW, 


—Ur 


-dw a 
u 


-dp a 

-7 p^r 
ox 

7 P d ( p a r ) 

r dr 

7 P dp a 


dx J ' r dcj) 
1 d (v a r) 1 dwa du a ' 

r dr r dcf) dx 


0 

0 

0 



x 0 ) 


(31) 

(32) 

(33) 

(34) 


With the Green’s functions known, the general solution is given by the convolution of 
the Green’s functions, or adjoint solutions, with the source terms. That is, 


P> (x,£) = - 



{Vi • [u a (X!,x,w)g s (X! ,tl)\ 


—q s (xi,G) Vi • [u a (xi,x, oj)] } x exp [— iu ( t — G)] dcodtid^i (35) 
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where 


and 


Vi 


d d 1 d \ 
(9rr | ’ dr i ’ rq dtfii ) 


(36) 


U„ (Xi.Xjjw) = [tl„, tlo.MlJ (Xi,X„,Oj) 


(37) 


However, the divergence of the adjoint velocity is known from the equation for the adjoint 
pressure p a , given by Eqn. (34). Also, q s is zero at the source point for the adjoint solution, 
so that 


p'(xd) 



,-dpa 

lUJpa + U — 
OX i 


(xi, x, lo) q s (xi, t\) exp [—iuj (t, — ti)] dujdt\dx.\ 


(38) 


or, 


p'{x,t) = ~ 



d _ d 

— + u - — ) {p a exp [-icj {t - ti)]} q s (xi, H) dudtidxi (39) 

Ot\ OX i 


It is important to note that, at this stage of the analysis, the convective derivative acts on 
the adjoint solution. However, integration by parts yields, 


p'(x,t) = 



\p a (xi,x, u) exp [-iu (t - ti)]] {q s (x x , H)} dudhdx .1 (40) 

JJti 


where the convective derivative is given by 

D _ ( d _ d \ 

Dt 1 ~ [dt! + U dxj ( ) 

The spectral density for the intensity is then obtained by forming the autocorrelation for 
the pressure and taking it’s Fourier transform, as given in Eqn. (5). The spectral density is 
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then found to be given by 


s (x, U ) = HI • • ■ JIJ Pa (x l , x, Wi ) Pa (x 2 , X, W 2 ) ^ {q s ( X i , *1 ) } {ft (x 2 , t 2 ) } 

x exp [—« (cji + UJ2) t + + iw 2 t 2 ] S (w — u>2 ) dh)\du2dt\dt2dx.\dx.2 ( 42 ) 

Again, as in models based on the acoustic analogy, it is necessary to make some assumption 
about the correlation function for the source terms. Based on experimental measurements 
of the two-point cross correlation of the axial velocity fluctuations in a fixed reference frame, 
Tam and Auriault 1 assume that 

(m {q " <*>•**» w 2 {q ‘ (X2 ' fe)! ) = *r xp {-J| Sr)2 + " 2 + } (43) 

where, 

£ = X\ — X 2 , ri = yi-y 2 , ( = zi-z 2 , t = t\ — t 2 ( 44 ) 

It should be noted again, as in Model I, that this is a major assumption and it is only 
supported indirectly by experimental measurements. In addition, the turbulent statistics 
are modeled in terms of the cross correlation of the convective derivative of the source. This 
is shown below to be a crucial choice in the formulation. 

Following the sequence of integrations given by Tam and Auriault 1 it is straightforward 
to show that 

S(x '" ) = ^////// P “ (X1,X ’^ )P “ <X2 ' X '“ ) ^I 

“l# } exp {-|[ + ~ b 2 + f2 ] } («) 

Tam and Auriault 1 argue that the difference between p a (xi, x, — co) and p a (x 2 ,x, uj) is a 
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simple phase factor such that 


Pa (xi, X, -uj) ~ p a (x 2 , x, oj) exp -i— cos 0 (xi - x 2 ) (46) 

a oo 

0 is the polar angle measured from the downstream jet axis. This implies that the sources 
have a limited spatial extent. That is, the length scale i s is small compared to the total 
extent of the source region. Thus, this model would not be appropriate for large scale or 
non compact sources. They also note that, 

Pa (x 2 , X, bj) = p* a (x 2 , X, Cd) (47) 


The integrations with respect to xi may now be replaced by integrations with respect to 
the separation distance £ = xi — x 2 The spectral density is then found to be given by, 


S (x, u) 



_ cxp{ .j-TTjPri) 

, S [1 + UJ 2 Tg (1 — U COS 9 /(loo)] 2 


In order to obtain a closed form result it is necessary to determine \p a (x 2 , x, uj)\ 2 . Since at 
9 = 7t/2 the effects of the mean flow are negligible, the equations for the adjoint functions 
(31) to (34) may be reduced to a Helmholtz equation for the adjoint pressure, 


V 2 p a + k 2 p a = x— -^<5 (x — x 2 ) (49) 

2ttc 2 0 

So that, 

9 C(J 2 

|P * (X2 ' X '"° I (50) 
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Then, from equation (48), with 9 = 7r/2, we obtain 


S (x, u) 


167 T 2 p 0 C%X 2 



exp{— ci o 2 l 2 / (4u 2 )} 
(1 + w 2 r 2 ) 


dx 2 


(51) 


With the scales given by Eqn. (23) and q 2 / c 2 = 4A 2 p 2 s k 2 /9, the contribution to the spectral 
density from an elemental volume of the turbulence is given by, 


dS (x, lo ) 


A A2 4 2 ,,7/2 (^’M) 2 

9p 0 CgX 2 4 fs ' (1 + 4 tt 2 uj 2 /uj 2 s ) 



7 r 2 c 2 w 2 A: 
c 2 cj 2 u 2 


dV 


(52) 


This is very nearly the same result as given by the Lighthill Acoustic Analogy in equation 
(24) except for the frequency weighting function. 


(Uj/idsf 

(1 + 47 t 2 uj 2 /uj 2 ) 



7 r 2 c 2 u) 2 k 1 

C 2 UJ 2 u 2 J 


(53) 


Figure 3 shows a comparison between predictions based on Tam and Auriault’s model and 
experiments. The agreement between the predictions and measurement is excellent. This 
excellent agreement follows that shown by Tam and Auriault 1 for a very wide set of jet 
operating conditions. The coefficients used in the present prediction based on Eqn. (52) are 
e T = 0.308, ci = 0.130, and A = 0.773. These are similar to the coefficients determined by 
Tam and Auriault 1 . Also shown in Fig. 3 is a prediction of the spectral density based on the 
coefficients determined by Tam and Auriault 1 . It can be seen that their predicted spectral 
density falls more rapidly at higher frequencies than the present predicted spectrum or the 
measurements by Tanna et al. 11 However, Tam and Auriault’s predictions agree well with 
measurements in other facilities. This shows how the values of the coefficients are influenced 
by the choice of experimental data. However, for both sets of coefficients, the values are 
arguably more reasonable than the value of c T obtained in Model I that was greater than 
unity. 
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4 Lighthill’s Acoustic Analogy: Model II 

In this section we will follow the same general approach used by Tam and Auriault: 1 but, 
here we apply the method to Lighthill’s acoustic analogy equation. The starting point is 
Eqn. (1). The Green’s function is given by the solution of the equation 


d 2 G 2 d 2 G 

dt 2 0 dxidxj 


8 (x — x x ) 8 (t — ti) 


(54) 


We seek a periodic Green’s function such that 


G (x, x-l, t, t\) = / G (x, x x , uj) exp [— iui (t — ti)] cluj 


(55) 


Then, the Fourier transform of equation (54) gives 


(V 2 + k 2 ) G (x, x x , u) = S ( 56 ) 

The adjoint function G a (x, x s ,uj) now satisfies the equation, 

(V 2 + A: 2 ) G a (x,x s ,u) = - S ^ 2jrc ^ (57) 


and it is readily shown that, 


G(x„x,,u j) = G, a (x^Xj.w) 


(58) 


Thus the far field density may be written 


f> (x,t) = 



G a (x 1? x,o ;) 


d 2 T ;■ 




dxudxij 


(x!, tj exp [— ico (■ t — G)] clujcltidhti (59) 
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The subsequent analysis is simplified if the derivatives are transferred onto the adjoint func- 
tion by integration by parts. That is, 


p' (*,£) 


d 2 G a 

dxudxij 


(x!, x, Lo)Tij (x,.. t x ) exp [—iuj ( t — ti)] dcodtidxi 


(60) 


The spectral density for the intensity is defined by Eqn. (5) and is given by, 


S (x,oj) = — 

Po 




d 2 G n 


( , d 2 G a 

(X"! , X, UJ \ ) 


dxudxij dx 2 k dx 2 1 

x exp [—iu> i ( t — t\ ) — ico 2 ( t — t 2 )] 5 (ui — ui 2 ) du]iduj 2 dt\dt 2 clx.idx. 2 


(x 2 , x, lo 2 ) (Tij (xi, t x ) T k( (x 2 ,t 2 )) 

(61) 


It is now assumed that the two point cross correlation of the Lighthill stress tensor in a fixed 
frame of reference is given by, 


(Tij (xi,t JTu (x 2 , t 2 )) 


A ijU p 2 u * exp 



i[K-,lr) 2 + y + C 2 ]} 


(62) 


Then, 


S (x,cj) = — 
Po 




d 2 G a , x d 2 G a , 

(Xi, X, Wi) — (x 2 , X, u 2 ) 


xA ijk£ p^4 exp 
+ {Vi ~ 2 / 2 ) 2 + (a ~ z 2 ) 2 ]} 


dxudxij 
\xi - x 2 


dx 2 k dx 2£ 


1 _ 2 
- yy ((aq - x 2 ) ~ u [h - t 2 )) 
ut, a 


x exp [— i (lj i + uj 2 ) t + + iuj 2 t 2 ] S (lj — ui 2 ) du)idu) 2 dtidt 2 dx.id'x. 2 (63) 


Again, following the sequence of integrations given by Tam and Auriault, 1 we obtain, 


S (x,cj) = 



Po 

( u 2 G s 
X6XP “TA 


d 2 G n 


dxudxij 


B 2 G f 

(xi, x, -u) - (x 2 , x,cj) A, ikl/ j 2 /P — 

dx 2 k dx 2l u 


ex p \ - i (?/ 2 + C 2 ) \ d*idx 2 


UTq 


u 


(64) 
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The solution to Eqn. (57) is 


G a (xi,x, -u) = 


1 exp [ik |xi — x| 
87 r 2 c 2 |x x — x| 


(65) 


It may be shown that as |xi — x| — ?• 00, 


dGg ^ f s dG a 

dxudxij ’ ’ dx 2 kdx 2 i 


(x 2 , X, OJ ) 


u 

647 t a c 3 x 2 


PiPjPkPt exp 


iu> (x 2 



( 66 ) 


where [l, = (xu —x,{)/x, etc. are direction cosines. Then, at 6 = 7 t /2 where, (x 2 — x x ) -x/rc ~ 
0, the far field spectral density for the intensity may be written, 


S (x,cj) = 


7T OJ 



,4 2 4 4 

^4 2222 p s u s — 


32n 3 p 0 c b x 2 

x exp (“1^1 exp I “Ir ~^r~^ ( r > 2 + C 2 ) rf * rfX2 


UT q 


U 


P, 


(67) 


With A 2222 = A 2 , and following the integrations with respect to £, we obtain 


S (x,cj) = 


7T 


,4 2 


167T 2 p 0 C b X^ 



faW i^Ts) 4 


Tf 


[1 


w 2 r 2 ; 


exp 


4fi 2 


(7x<; 


( 68 ) 


Using the scales defined by Eqn. (23) and with u s = y/2k/3, the contribution to the spectral 
density from each volume element of the flow may be written, 


dS (x,cj) 


7T 2 y/7T A 2 G 2,7/2 (w/u; s ) 4 / 7T 2 C» CJ 2 A) \ 

7 — — 7 ~P 2 s k — , 9 9 , 9 . exp 5-^—^— dV 

9 p 0 c b x 2 c 3 (1 + 4tt 2 uj 2 /uj 2 ) \ c 2 4 x> 2 u 2 J 


(69) 


This is essentially the same result obtained from Tam and Auriault’s model 1 except for the 
slight but significant change in the frequency weighting factor 


( uj/lo s ) 4 / 7 r 2 e 2 oj 2 k \ 

(1 + 4tt 2 uj 2 /uj 2 ) ^ \ c 2 4c o 2 u 2 ) 


(70) 
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Figure 4 shows a comparison of predictions made with all three models with experimental 
data. The coefficients used in the Model II prediction are, c T = 1.0, q = 0.78, and A = 
0.8475. Clearly, the Model II prediction is an improvement over that made with the Model 
I version based on the acoustic analogy. Some of this is due to the additional freedom to 
specify three independent coefficients unlike the two available for Model I. However, the 
shape of the spectrum is still not as well predicted as that given by Tam and Auriault’s 
model 1 for the same level of empiricism. 

In the next section, we consider the final difference between the various models: the 
specification of the turbulent source statistics. 

5 Models for the Turbulent Source Statistics 

In the previous sections it was shown how a model based on the acoustic analogy failed 
to provide as good a set of noise predictions as the Tam and Auriault model 1 even when 
the problems were formulated in the same manner. However, there remains one issue to 
be considered: how the statistics of the turbulent sources are described. Tam and Auriault 
model the source statistics with Eqn. (43) whereas the model of the last section used Eqn. 
(62). In the former case it is assumed that the cross correlation of the convective derivative 
of the source term follows the experimentally measured cross correlation of the axial velocity 
fluctuations in the jet. In the latter case, it is assumed that it is the cross correlation of 
the Lighthill stress tensor that follows the measured form. In this section the two prediction 
schemes are compared when consistent forms of the statistics of the turbulent source are 
used. 

There are two ways to modify Tam and Auriault’s model 1 to make their assumptions 
consistent with Model II. The first is to retain the convective derivative operators on the 
adjoint solution in Eqn. (38). Once the form of the adjoint solution is obtained, the operators 
may be applied in a straightforward manner. However, in the general case, when the adjoint 
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solution is obtained numerically, it is as convenient to still perform the integration by parts 
and use the form for the spectral density given by Eqn. (42). However, it is now assumed 
that the cross correlation of the source function may be written, 


(Qs ( xi ,ti)q a (x 2 ,t 2 )) = 7! ex P<! 4 [it ~ ut ) 2 + rf + ( 2 ] 


(71) 


This is consistent with the assumption made in Eqn. (62) for Model II based on the acoustic 
analogy. The cross correlation required in Eqn. (42) is, 


W , (Xl ' t,] dT 2 (X2 ' h] ) = d^dT 2 <«> (x " tl) (X2 ' k)) 


(72) 


The derivatives may be written in terms of the time delay r and the axial spatial separation 
distance £ such that, 


D 2 


Dt\Dt2 


d d 

a7 + “^ 


(73) 


Now, if 


d d 

K - St) then ( I f = 0 


(74) 


Thus, 


Dq s , , Dq s 

x 1 ,t 1 )—-(x 2l t 2 ) 


Dti 


Dto 


Qs -2 

~ — 77 u exp 

c 2 


-J 2 [({ - urf + ,f + c 2 ] 


X 


d 2 


exp 


lei 


UTg 


(75) 


If this result is substituted into Eqn. (42) and the integrations with respect to co i, w 2 , H and 
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t-> are performed, we obtain, 


S (x, ui) 



d 2 


exp 




r/x | r/x 2 


( 76 ) 


The adjoint solutions may be replaced by the form given by Eqn. (50). Then the integrations 
with respect to xi may be replaced by integrations with respect to £. The integrations with 
respect to r) and £ are readily performed as before. However, the integration with respect to 
£ now takes the form, 



— OO 


This integral may be rewritten in terms of a new variable #'= £/ ( ur s ), such that, 


I 



,zj)] exp (— iuiT s z ) dz 


(78) 


Since the first term in the integrand must be treated as a generalized function, due to its 
discontinuous behavior at £ = 0, integration by parts should be performed before the integral 
may be evaluated. This is discussed in more detail in an Appendix to this paper. This gives, 


I 


d 2 

— / exp (- ^l) — [exp (-iur s z)] dz 


UT. 


uj 2 t s 


u 


exp (— |^| — iujT s z ) dz 


(79) 

(80) 


The integral may now be evaluated by separating the integral into the contribution from 
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oo < £ < 0 and that from 0 < £ < oo. This yields, 


2c o 2 r s 

u (1 + l o 2 r 2 ) 


( 81 ) 


The use of Eqn. (50) for the magnitude of the adjoint pressure then gives the final result for 
the spectral density, 


S (x, cj) 


Vtt 1 

167 r 2 poC^x 2 


c 2 r 3 (1 + cj 2 t 2 ) 


exp 


4 u 2 ) 


rlx 2 


(82) 


If g 2 /c 2 is set equal to M 2 p 2 ?4 this is exactly the result given by Eqn. (68) for Model II based 
on the acoustic analogy described in the previous section. 

Thus it may be stated that: both Model II based on the acoustic analogy and Tam and 
Auriault’s model 1 yield identical noise prediction formulas at 90 degrees to the jet axis if 
consistent assumptions are made in the statistical description of the turbuleiit sources. 


6 Discussion 

In the previous section it was shown that both a model based on an acoustic analogy and Tam 
and Auriault’s model 1 yield identical noise prediction formulas if consistent assumptions are 
made concerning the statistical description of the turbulent noise sources. This might appear 
to be a somewhat negative result, as it shows that the recent model by Tam and Auriault 1 
would give as poor a prediction as models based on the acoustic analogy for the same source 
description. However, it could be viewed as providing guidance on what assumptions should 
be made in the source description if good predictions are to be made. One could obtain 
a noise prediction formula identical to that given by Tam and Auriault 1 in Eqn. (52) in a 
number of ways. For example, only two of the spatial derivatives acting on the Lighthill 
source term correlation in Eqn. (59) need to be transferred to the Green’s functions. Then 
an assumption could be made that the two-point cross correlation of the new “source” term 
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is given by 


dTjj dTk( 


dXi 


dxi 


(x 2 ,t 2 ) ) = V/A exp 




c 2 ] (83) 


However, this is just as arbitrary an assumption as that made by Tam and Auriault 1 in 
proposing the form of two-point cross correlation given by Eqn. (43). Alternatively, we need 
to find a form of two-point cross correlation function for the Lighthill stress tensor that would 
give the form assumed by Tam and Auriault 1 when the convective derivative operations are 
applied, as in Eqn. (72). This is no easy matter because of the discontinuous derivative of 
the assumed correlation at £ = 0. However, as noted by one reviewer, |£| could be replaced 
by a /£ 2 + S 2 and letting 8 — >■ 0 in the final result. 

In the remainder of this final section we propose an acoustic analogy formulation that 
includes, as a part, the source description proposed by Tam and Auriault 1 : but, it also 
identifies other physical mechanisms that are equally good candidates for noise sources. 
Other formulations of this type could be proposed: however, our formulation is based on 
equations of motion written in terms of the logarithm of the pressure. 

The equations of continuity, momentum and energy, and the equation of state for a 
perfect gas, may be rearranged in the form, 


Dn dui 

Dt Ox,, 


and 


Du, 2 dn 
Dt, dx, 

where 



(84) 


(85) 


( 86 ) 
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( 87 ) 


d _ a a 

Dt dt ' O.i'j 

and p 0 is the mean static pressure that is assumed to constant. Equations (84) and (85) may 
be recast in the form of a third-order equation [Lilley , 9 Colonius et a/. 15 ], 


D 

'D 2 7T 

d 

(c 2 —\ 

„ dm 

' d 

A— y 

Dt 

Dt 2 

dx i 

V Sx,)\ 

dxj 

.dxj 

V 8xJ\ 


0 duj du k duj 
d Xj dx.j dx k 


( 88 ) 


If the variables in this equation are linearized about the mean thermodynamic properties 
and a parallel mean flow, given by, 


u (x, t) = U {x 2 , x 3 ) S n + u'i (x, t) 


(89) 


Lilley’s equation 9 may be obtained in the form 


Do 

Dt 


Dt 2 



+ 2 c 2 


dU d 2 7T' 
dx 2 dxidx 2 


dU dV 
dx 3 dx | dx 3 


r 


(90) 


A detailed expansion of T is given by Colonius et a /. 15 All terms in T are at least second 
order in the fluctuations. In the limit of infinitesimal disturbances the equation reduces to 
a homogeneous equation that describes the propagation of sound in a parallel shear flow, 
such as that developed by Pridmore-Brown . 16 Also, as noted by Colonius et a /., 15 it is in the 
form of an acoustic analogy since it is equivalent to the equation that describes the pressure 
fluctuations generated by an external distribution of stresses (and other source terms that 
are often neglected) imposed on a parallel shear flow. 

Alternatively an acoustic analogy could be equally well formed by simply splitting the 
variables in Eqns. (84) and (85) into their mean and fluctuating components, assuming a 
parallel mean flow, and retaining only linear terms on the left hand side of the equation. 
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This yields, 


Dy 

Dt 


d< 

dxi 


= — w 


( 91 ) 


and 


with 


DqK 

Dt 


( dU 
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\dx' 


~Un 


clU 

dx? 


- u i 
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T dn' 


i 1 


2 _____ _ 


( 9 : 7 ’, 


( 92 ) 


A, 

Dt, 



d 

dx\ 


( 93 ) 


and 


// / 

» = 

J OX ; 


and 


U ' Ox ; 


2 ,(hd_ 

dx; 


( 94 ) 


The double prime notation is used to denote terms that are second order in the fluctuations. 
The term w" represents the net rate of work done at second order by the fluctuating pressure 
on a fluid element. The term /" is the same as that obtained by Goldstein, 1 ' being equivalent 
to an externally applied force per unit mass acting on a parallel sheared mean flow. This 
term may also be written, 


= — £. 


i jkUjU) k 


d_ 

dx; 




2 ' 

dx,; 


( 95 ) 


where is the alternating tensor and ui' k is the vorticity tensor. The first term in the total 
fluctuating force (per unit mass) is called the “vortex force” (also written in vector notation 
as u' x u>'). It is associated with the transport of rotating fluid elements by the fluctuating 
velocity field. The second term is the force generated by the second order dynamic head. It 
should be noted that this term is essentially the same single source term proposed by Tam and 
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Auriault . 1 The last term in Eqn. (95) is a source associated with temperature fluctuations in 
the turbulent flow as discussed by Tester and Morfey . 18 Note that c 2 ' = 7 RT' . This source 
term is an important component of the semi-empirical noise prediction model developed by 
Morfey et a /. 19 This last term appears as a source in the framework of an acoustic analogy. 
However, it is readily identified as the effect of sound propagation through a medium with a 
fluctuating speed of sound. 

Equations (91) and (92) could be rearranged into a single inhomogeneous equation for 
7 r', exactly in the form of Eqn. (90). However, in their unmanipulated form it is easier to see 
the equivalent sources in the momentum and energy equations. The source term F, given by 
Colonius et ah , 15 is reproduced exactly by this rearrangement 1 . (To obtain the form given 
by Colonius et a /., 15 the source associated with speed of sound fluctuations is neglected as 
are third and higher order terms.) The complexity of F is formidable. In its simplified form, 

r w d q f du >j \ o f dU duty ] dU duty \ (gg) 

Dt \dxidxj J \clx2 dx\dxj <lx$ dx\dxj J 


This has been interpreted as a quadrupole source distribution with components referred to 
as self- and shear-noise. However, the form of source term given by Eqn. (96) is simply 
the result of recasting the original inhomogeneous, linearized Euler equations given by Eqns. 
(91) and (92). In that primitive form there is no source term, or by inference, no source 
mechanism, associated with the mean velocity gradient, or any other mean flow property. All 
equivalent source terms are second order in the fluctuations only. In summary, the acoustic 
analogy formed by Eqns. (91) and (92) describes acoustic radiation from a sheared mean 
flow due to externally-applied forces and work performed on the fluid. 

Since the source terms on the right hand sides of Eqns. (91) and (92) are not all in 


1 There appears to be a minor typographical error in the final source term component Vb, that should be 


On 1 ( , du'2 \ 
dx\ \ * dxj J 


+ U-2 


d 

dxi 
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the form of the divergence of a scalar, as in Tam and Auriault’s model, 1 it is not possible 
to follow the simplifying analysis of Eqns. (35) to (40). In addition, since there is more 
than one source term, their correlation must be considered when the autocorrelation of the 
far field pressure is formed. Previous analyses have considered each source mechanism to 
be statistically independent and that approach, though somewhat questionable, could be 
followed. 

The scaling of the source terms, based on the CFD solution, may also be addressed. The 
first two components of /" of Eqn. (95) should be scaled by 

k/t a (97) 


The last term should be scaled by 


- T 0 ) p s k 

T 0 p 0 4 

where, T s is the mean temperature in the source volume element. It should be noted that this 
is only distinguished from the scaling for the first source terms (or the source term in Tam 
and Auriault’s model 1 ) by the local relative temperature difference. Finally, the unsteady 
work source term in Eqn. (91) should scale with 

1 p s A: 3 / 2 
c 2 po 4 

There are several issues that remain to be addressed. Though the present paper has 
considered the noise radiation at 90° to the jet axis, so that mean flow/acoustic interaction 
effects may be neglected, there are distinct differences between the traditional (Model I) 
formulation and Tam and Auriault’s model 1 at other angles to the jet axis. In the latter 
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model, the variation of spectral density with polar angle is contained in the factor 

\Pa{x 2 ,X,u)\ 2 

[l + (w 2 / w s) (1 ~ ucos9/c 0 ) 2 ] 

Away from the zone of silence there is little angular variation in \p a (x2,x,uj)\ at a given 
frequency as shown by Tam and Auriault. 14 The factor in the denominator in Eqn. (100) is 
a very weak function of polar angle. For example, with u/c 0 = 1.0 there is only a 2 decibel 
change in the spectral density for uj / oj s = 1.0 and 6 = 60° relative to the value at 90°. This 
is consistent with the experimental spectra shown by Tam and Auriault. 14 However, the 
traditional models based on Lighthill’s acoustic analogy suggest a directivity that varies as 
(1 — M c cosd)~ 5 . If we let M c = 0.62 Uj/c 0 (Goldstein' ) then a relative change of 8 decibels 
is predicted for the same conditions. This is much greater than the measured variation. Fur- 
ther analysis by Goldstein' to include the mean flow effects changes the convection factor 
from (1 — M c cosd)~ 5 to (1 — M c cos#) -3 . This results in predictions for the directivity of 
the overall intensity in better agreement with experiment. In addition, in traditional models 
based on the acoustic analogy, there is no separation of source mechanisms between con- 
tributions from small-scale or large-scale turbulence. So predictions are made at all polar 
angles for the same assumed source. Whereas, in Tam and Auriault’s model, 14 predictions 
are limited to angles approximately greater than 60° to the jet downstream axis. Closer to 
the jet axis Tam and Auriault argue that a noise mechanism associated with the large-scale 
turbulence is present. It should be noted that Model II based on the acoustic analogy, de- 
scribed in Section 4, gives a convection factor identical to that of Tam and Auriault. The 
reason for the difference is tied to the choice of two-point cross correlation factor for the 
sources. In Model I the form of the cross correlation is assumed in a moving frame of ref- 
erence. Whereas, in Model II, the form of the cross correlation is assumed in a reference 
frame fixed to the jet nozzle. The fact that the two models give different convection fac- 
tors indicates that one or other of the assumed forms is incorrect or insufficiently accurate. 


30 



The moving and fixed reference frame cross correlations are related simply by a coordinate 
transformation. So consistent choices are possible and the final result must be independent 
of the reference frame. It is readily shown that the usual Gaussian form of cross correlation 
in the moving frame gives a fixed frame correlation that is consistent with Eqn. (71) and 
measurements. Since, as noted above, the temporal variations of the sources that determine 
the sound radiation are more clearly identifiable in a moving reference frame, it appears to 
be the appropriate choice: and with its selection comes convective amplification. Finally it 
should be noted that the measured two-point correlations should contain contributions from 
both the large-scale and small-scale turbulence. The wavenumber /frequency spectrum that 
provides the distinction between radiating and non-radiating components of the turbulent 
field [see Eqn. (17)] should be able to describe both components if it is chosen appropriately. 
The extraction of two-point correlations from either experimental or simulation data for 
high speed jets, where it known that the instability wave radiation is present, would be very 
illuminating. Thus, the issues of the appropriate convection factor and the separate contri- 
butions, if they are separate, of the small-scale and large-scale noise generation mechanisms 
remain open questions. 

In this paper we have tried to reconcile the apparent differences between models based 
on the acoustic analogy and alternative models. We have shown that there is no difference 
in the eventual noise prediction formulas if consistent assumptions are made concerning the 
statistical properties of the turbulent noise sources. We have also proposed an acoustic 
analogy that provides a clear representation of equivalent source mechanisms. We have 
also tried to identify some of the remaining issues that have yet to be resolved and have 
proposed measurements or simulations that would help to resolve these issues. The recent 
reexamination of models based on the acoustic analogy, much of it prompted by the success 
of the model proposed by Tam and Auriault 1 has been very valuable. The authors hope that 
this paper will stimulate further work that will result in an improved understanding of how 
to predict radiated noise using CFD simulations. 
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7 Appendix. The Use of Generalized Differentiation 
in the Algebraic Manipulation of Equation (78) 


Consider the following integral depending on the variable y and a real variable ft 


OO 

/(,) = /“ p( -'*- v l+mdx 
— OO 

2 exp (ijjy) 

1 + 


( 101 ) 

( 102 ) 


It is clear that the second derivative of the function on the left hand side of Eqn. (102) is 


d 2 f ^ = 2ft 2 exp jifty) 


dy 2 


1 + ft 2 


(103) 


If we differentiate Eqn. (101) with respect to the variable y we obtain 

OO 

r d 2 

/ -r-z {exp (— \x — y I + //3.x)} dx 
J dy 2 

— OO 

It might appear to be possible to split this integral into two parts giving, 

y oo 

f d 2 , , . . , f d 2 , , . ,, , 2exp(ifty) 

/ “i (7 {exp (.x - y + iftx)\dx + / — - {exp (-.x + y + %ftx)\ dx = — — — (lOo) 

J dy z J dy z 1 + ft z 

-oo y 

This does not agree with the result given by Eqn. (103). To determine why this is we 
note that the second derivative of the function exp (— \x — y |) is discontinuous at x = y and 
thus taking the second derivative inside the integral is not valid. We should have kept the 
derivatives outside the integral, split the range of integration, and then used the Leibniz rule 



32 



of differentiation under the integral sign twice to get, 


d 2 

dy 2 


exp (— \x — y\ + //3a’) da 


2 ^ / exp (a — y + ?'/3a) dx + j exp (—x + y + i,6x) dx 


2/3 2 exp ( i >'//) 
1 + /3 2 


(106) 


We have obtained the correct result but only at the expense of more algebraic manipulations. 
We want to be able to obtain the correct result when we change the order of differentiation 
and integration without the added burden of book-keeping when our integrand is discontin- 
uous. This can be achieved as follows. 

It is seen that f (y) is an analytic function. This means that ordinary and generalized 
derivatives of this function are identical. However, we can exchange the order of integration 
and generalized differentiation. Using a bar over the differentiation sign to denote generalized 
differentiation, we have now the following mathematically correct manipulation: 
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dy 2 


r d 2 

( y ) = J j ~2 { ex P (- l x - V I + ifa)} dx 


— oo 
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{exp — | a — y |) e lBx dx 


— OO 
OO 


I exp — \x — y\ -j-r; {e t(3x } dx 

o 

2/3 2 exp ( i in) 

TTp 2 


(107) 


So we have obtained the same result as Eqn. (103). 

In this paper we have obtained an equation similar to Eqn. (104) and then we claimed 
that the second derivative with respect to the variable y must be treated as a generalized 
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derivative. The purpose of this appendix is to justify this claim. The root of the problem in 
Eqn. (104) is that the interchange of the order of limit operations of integration and ordinary 
differentiation is not permissible because of the discontinuity of the second derivative of the 
integrand at x = y. However, in this paper, we have freely interchanged the order of 
integration and ordinary differentiation without any attempt at justification. For example, 
consider the derivation of the relation, 


d 2 T xx 

dt 2 

d 2 T \ 

(yi,* i} -^f (y 

d A 

= dr 4 i 1 (yii t) T XX ( y 2 • To)) 

(108) 


To 

x (y 2 — Vi) 
= t + T+~- — — 

X c 0 

(109) 


This is necessary in proceeding from Eqn. (6) to Eqn. (9) in the paper. We have performed 
two integrations by parts with respect to t, replaced the fourth order differentiation with 
respect to t with the same order differentiation with respect to r, and then we exchanged the 
order of differentiation with the integration for time averaging (see Goldstein, 7 Appendix 2). 
These operations can only be justified for integrands whose derivatives are continuous. More 
precisely, the uniform convergence of the integrals must be considered. As has happened 
in this paper, we may next introduce an algebraic expression for the statistical properties 
of the turbulence and perform algebraic manipulations that involve bringing into integrals 
derivatives of some order of a variable. We can do this easily and without the burden of 
the detailed analysis of classical calculus, exhibited in the manipulations of Eqn. (106), if 
we treat all derivatives from the start as generalized derivatives. Practically all the analysis 
remains the same except in situations like the exchange of the order of limit operations 
where we must indicate generalized differentiation explicitly as we have done here. The 
situation is completely analogous to linear aerodynamic theory where occasionally we have 
to interpret divergent integrals as the finite part of divergent integrals. The reason is that 
divergent integrals appear if improper exchange of the order of differentiation and integration 
is performed. The finite part procedure gives the same result as if the derivative was kept 
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outside the integral. The procedure itself is obtained by assuming that the derivative is a 
generalized derivative and can, therefore, be taken inside the integral and interpreted by the 
rules of generalized function theory. 

Readers should consult Canonical Regularization in Gelfand and Shilov , 20 for the theo- 
retical background relevant to the subject of this appendix. 
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Figure 1: Comparison With Experimental Data of the Predicted Spectral Density at 90 
to the Jet Axis Using Model I Based on the Acoustic Analogy, x = 72 Dj, Mj = 0.911 
Tj/T 0 = 0.975, Dj = 0.0508 m. Q O • Prediction; • , Experimental, Taima et al. 11 




Figure 2: Prediction of Spectral Density at 90° to the Jet Axis Using Model I Based on the 
Acoustic Analogy Showing the Contribution to the Spectrum in One Diameter Slices. The 
Contributions From Slices for the First Twenty Diameters are Shown in Alternating Solid 

and Dashed Lines, x = 72 Dj, Mj = 0.911, Tj/T 0 = 0.975, Dj = 0.0508 m. O Cb 

Overall Prediction. 
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Figure 3: Comparison With Experimental Data of the Predicted Spectral Density at 90° to 
the Jet Axis Using Tam and Auriault’s Model. 1 x = 72 Dj, Mj = 0.911, Tj/T 0 = 0.975, 

Dj = 0.0508 m. □ □ , Present prediction using coefficients determined by comparison 

with Tanna et al. 11 data; O O, Prediction based on coefficients determined by Tam 

and Auriault; 1 • , Experiment, Tanna et al. 11 
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Figure 4: Comparison With Experimental Data of the Predicted Spectral Density at 90° to 
the Jet Axis Using Models I and II Based on the Acoustic Analogy and Tam and Auriault’s 
Model. 1 x = 7-D j, Mj = 0.911, Tj /T 0 = 0.975, D, = 0.0508 m. O O, Model I; 

A A, Model II; □ □, Tam and Auriault’s Model; • , Experiment, Tanna et al. 
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